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ABSTRACT 

Based on kinematic data observed for very young, massive clusters that appear to be in dynamical 
equilibrium, it has recently been argued that such young systems set examples where the early residual 
gas-expulsion did not happen or had no dynamical effect. The intriguing scenario of a star cluster 
forming through a single starburst has thereby been challenged. Choosing the case of the R136 
cluster of the Large Magellanic Cloud, the most cited one in this context, we perform direct N-body 
computations that mimic the early evolution of this cluster including the gas-removal phase (on a 
thermal timescale). Our calculations show that under plausible initial conditions as consistent from 
observational data, a large fraction (> 60%) of a gas-expelled, expanding R136-like cluster is bound 
to regain dynamical equilibrium by its current age. Therefore, the recent measurements of velocity 
dispersion in the inner regions of R136, that indicate that the cluster is in dynamical equilibrium, are 
consistent with an earlier substantial gas expulsion of R136 followed by a rapid re-virialization (in w 1 
Myr). Additionally, we find that the less massive Galactic NGC 3603 Young Cluster (NYC), with a 
substantially longer re-virialization time, is likely to be found deviated from dynamical equilibrium at 
its present age (w 1 Myr) . The recently obtained stellar proper motions in the central part of the NYC 
indeed suggest this and are consistent with the computed models. This work significantly extends 
previous models of the Orion Nebula Cluster which already demonstrated that the re-virialization 
time of young post-gas-expulsion clusters decreases with increasing pre-expulsion density. 
Keywords: stars: kinematics and dynamics — open clusters and associations: individ- 
ual(R136,NGC3603) — galaxies: star clusters: general — galaxies: individual(LMC) 



; 1. INTRODUCTION 

The question of how a young star cluster is formed 
has been debated for decades. According to the classi- 
cal scenario, a cluster forms essentially through a sin- 
gle star-burst event. The individual proto-stellar cores 
within the parent or proto-cluster gas cloud approach 
their hydrogen-burning main-sequences (MS) to form an 
infant star cluster which still remains embedded within 
the residual gas that did not collapse to stars. This resid- 
ual gas receives kinetic energy and radiation-pressure 
from the radiation of the massive MS and pre-main- 
sequence (PMS) stars until it gets unbound from the sys- 
tem and then escapes. This gas-removal process can be 
expected to be very rapid — typically faster than or sim- 
ilar to the dynamical crossing time of the embedded clus- 
ter (Lada & Lada 2003). The remaining gas- free cluster 
must expand due to the corresponding dilution of grav- 
itational potential well. For the hypothetical case of in- 
stantaneous gas removal, the resultant cluster should get 
unbound if the total mass lost as gas is equal to or more 
than the mass remaining in the stars (i.e., star formation 
efficiency (SEE) e < 50%), as is true for any gravitation- 
ally self-bound system. For slower gas removal, the sur- 
vivability of the gas-deprived cluster as a bound system, 
for e < 50%, increases (Kroupa 2008). Two-body relax- 
ation, which evolves the cluster towards a higher central 
concentration until the beginning of gas-expulsion, also 
enhances survival. 

The above scenario, applicable to the formation of 
Galactic and extra-galactic young star clusters, has been 
first depicted by, e.g., Hills (1980), Elmegreen (1983) 
and Mathieu (1983) either analytically or by small N 



direct N-body studies (Lada et al. 1984). These mile- 
stone works have later been elaborated by Kroupa et 
al. (2001), Geyer & Burkert (2001) and Baumgardt & 
Kroupa (2007) by extensive direct N-body calculations 
of gas-expelling model infant clusters. Such theoretical 
studies, along with observations of young embedded sys- 
tems (Lada 1999; Elmegreen et al. 2000), suggest a min- 
imal star- formation efficiency of e « 1/3 for forming a 
bound, gas-free cluster, applicable to embedded systems 
initially in dynamical (or virial) equilibrium. An observa- 
tional evidence of the residual gas expulsion process can 
be the < 2.5 Myr old Orion Nebula Cluster (ONC) whose 
measured velocity dispersion implies that it must be ex- 
panding (Jones & Walker 1988). While the ONC was 
widely considered to end-up in an unbound OB associa- 
tion {e.g., Zinnecker et al. 1993), the detailed modelling 
by Kroupa et al. (2001) implies that it is rather destined 
to become a bound cluster very similar to the Pleiades 
by 100 Myr. Recent proper motion measurements of the 
Galactic w 1 Myr old NGC 3603 Young Cluster (NYC) 
by Rochau et al. (2010) also suggests that its stars might 
be away from energy-equipartition. Notably, a demon- 
stration of the expansion of clusters younger than a few 
Myr as they age is collated by Brandner (2008). 

The key properties of the stellar initial mass function 
(IMF; Kroupa et al. 2013), in particular, its observed uni- 
versality is also supported by a monolithic cluster forma- 
tion scenario that involves competitive gas accretion by 
the most massive proto-stellar cores. For example, de- 
tailed three-dimensional, adaptive-mesh FLASH calcula- 
tions including radiation feedback by Peters et al. (2010, 
2011) well reproduce not only the general form of the uni- 
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versal IMF but also the relation between the total mass 
(in stars) of an cmbcckled cluster and the mass of its 
most massive stellar member (Weidner & Kroupa 2004). 
Recently, Liu et al. (2012) find observational signatures 
of competitive accretion in the proto-cluster GlO.6-0.4. 

Very recently however, questions have been raised 
against the clustered mode of star formation, preferring 
instead a hierarchical or continuous star formation {e.g., 
Bressert et al. 2010; Gutcrmuth et al. 2011). The lat- 
ter is inferred from the stellar density distribution in 
young, star forming regions over a "global" scale. No- 
tably, Pfalzncr et al. (2012) showed that unless the den- 
sity contrast among the individual clusters is very high 
any such conclusions based on surface density profiles 
can be highly ambiguous and is useful only when suffi- 
ciently deep and complete data are available. Further- 
more, from the kinematics of several very young clus- 
ters indicating their dynamical equilibrium, in particular 
that of the R136 cluster in the Large Magellanic Cloud 
(LMC), it has been argued that they must have avoided 
any substantial residual gas-expulsion phase {i.e., effec- 
tively have had 100% star-formation efficiency, also see 
Goodwin 2009 in this context); otherwise they would 
have been found expanding at such a young age. This has 
been particularly emphasized very recently by Henault- 
Brunct et al. (2012). A key ingredient that might be 
missing in such arguments is the consideration of the re- 
virialization of the expanding gas-expelled cluster, i.e., 
quick formation of a bound system after a fraction of the 
expanding cluster reverses back on a free-fall timescale. 

The possible fate of the expanding ONC being a 
Pleiades-like bound system has been shown by Kroupa et 
al. (2001) whose models already demonstrate that the re- 
virialization time decreases with increasing initial cluster 
density (Fig. 1 of Kroupa et al. 2001). The survivability 
as a bound cluster after gas-removal has been studied 
later in detail by, e.g., Baumgardt & Kroupa (2007). In 
this work, we find it crucial to study this issue again 
with particular focus on the case of R136. This becomes 
essential because of the current interpretations of the 
kinematic data of young, massive clusters such as R136, 
that seem to contradict the classical picture of monolithic 
cluster formation and the associated event of residual 
gas expulsion. We show that under plausible conditions, 
the re-virialization of an R136-like massive cluster, after 
the gas-removal, is prompt enough to make it be found 
in virial equilibrium despite of its young age. The re- 
cently obtained low (line-of-sight) velocity dispersion of 
single 0-stars in R136 (Henault-Brunet et al. 2012) is 
also found to be consistent with a re-virializcd system. 
We also find that the much lighter NYC, on the other 
hand, is unlikely to be virialized at its current age of ~ 1 
Myr (Stolte et al. 2004, 2006), as indicated by Rochau 
et al. (2010). 

2. MODEL COMPUTATIONS 

2.1. Gas removal 

A thorough modelling of gas-removal from embedded 
clusters is complicated by the radiation hydrodynamical 
processes which is extremely complex and involves un- 
certain physical mechanisms. For simplicity, we therefore 
mimic the essential dynamical effects of the gas-expulsion 
process by applying a diluting, spherically-symmetric ex- 



ternal gravitational potential to a model cluster as in 
Kroupa et al. (2001). Specifically, we use the potential 
of the spherically-symmetric, time(f)-varying mass dis- 
tribution 

Mg{t)= Mg{0) t<Td, 

Mg(i)=M,(0)exp(-(^^)i>Td. (1) 

Here, Mg{t) is the total mass in gas which is spatially 
distributed with the same initial Plummcr density distri- 
bution (Kroupa 2008; sec below) as the stars and starts 
depleting with timescale Tg after a delay of r^. The Plum- 
mer radius of the gas distribution is kept time-invariant 
(Kroupa et al. 2001). Such an analytic approach is par- 
tially justified by Geyer & Burkert (2001) who perform 
comparison computations treating the gas with the SPH 
method. 

The exact values of the essential parameters quantify- 
ing the gas-expulsion timescale, viz., Tg and depend 
on gas-physics. For simplicity, we use an average gas ve- 
locity of Vg K, 10 km s~^ which is the typical sound-speed 
in an H II (ionized hydrogen) region. This gives 
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where r/i(0) is the initial half-mass radius of the stel- 
lar cluster/gas. The coupling of stellar radiation with 
the ionized residual proto-cluster gas over-pressures the 
latter substantially and can even make it radiation- 
pressure-dominatcd (RPD), for massive clusters in the 
initial phase of the expansion of the gas. During the 
RPD phase, the gas is driven out at speeds consider- 
ably exceeding the sound-speed of the ionized medium 
(Krumholz & Matzner 2009). Once the expanding gas 
becomes gas-pressure-dominated (GPD), it then contin- 
ues to flow out with the sound-speed of an H II region 
(Hills 1980). Hence, the above Tg, from Vg w 10 km s^^, 
represents its upper limit; it can be shorter depending 
on the duration of the RPD state (also see Sec. 4). No- 
tably, the initial RPD phase is crucial to launch the gas 
from very massive systems whose escape speed exceeds 
the sound-speed in H II gas (Krumholz & Matzner 2009). 

As for the delay-time, we take the representative value 
of Td ~ 0.6 Myr (Kroupa et al. 2001). The correct value 
of Td is again complicated by radiative gas-physics. An 
idea of Td can be obtained from the lifetimes of the ultra- 
compact H II (UCHII) regions which can be upto w 10'"' 
yr (0.1 Myr; Churchwell 2002). The very compact prc- 
gas-expulsion clusters (Sec. 2.2) have sizes {rh{0)) only a 
factor of f» 3 — 4 larger than the typical size of a UCHII 
region (w 0.1 pc). If one applies a similar Stromgren 
sphere expansion scenario (Churchwell 2002 and refer- 
ences therein) to the compact embedded cluster, the es- 
timated delay-time, Td, before a sphere of radius r^(0) 
becomes ionized, would also be larger by a similar fac- 
tor and hence close to the above representative value. 
Once the gas is ionized, it couples strongly with the ra- 
diation from the 0/B stars and launched immediately 
(see above) . High- velocity jet outflows from proto-stars 
(Patel et al. 2005) additionally facilitate the gas removal. 

For super-massive clusters (> IO^Mq), however, a 
"stagnation radius" can form within the embedded clus- 
ter inside which the radiation cooling becomes suffi- 
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Table 1 

Computed model parameters. The ONC-A/B computations arc from Kroupa et al. (2001) and included for comparison. The quoted TvirS 
correspond to the re-virialization times of 30% of the initial cluster mass (in stars) after the delay time . The "BSE" column indicates 

the presence of stellar evolution. 



Cluster 


Mecl(O)/Af0 




rh{0)/pc 


Z/Zq 


T,/Myr 


Ter(0)/Myr 


Td/Myr 


BSE 


Tvir/Myr 


R136 


1.0 X 10^ 


2.0 X 10^ 


0.45 


0.5 


0.045 


0.021 


0.0, 0.6 


Yes 


0.9 


NYC 


1.3 X 10"^ 


2.6 X 10^ 


0.34 


1.0 


0.034 


0.038 


0.0,0.6 


Yes 


2.2 


ONC-A 


3.7 X 10=* 


7.4 X 10^ 


0.45 


1.0 


0.045 


0.23 


0.6 


Yes 


> 10 


ONC-B 


4.2 X 10^ 


8.4 X 10^ 


0.21 


1.0 


0.021 


0.066 


0.6 


Yes 


«3 



ciently efficient to possibly form second-generation stars 
(Wiinsch ct al. 2011). Also, as discussed above, the 
gas-outflow can initially be supersonic which generates 
shock-fronts. Although shocked, it is unlikely that star- 
formation will occur in such an RPD gas. Later, during 
the GPD outflow, the flow can still be supersonic in the 
rarer/colder outer parts of the embedded cluster where 
the average sound-speed might be lower than that typical 
for H II gas. However, it is not clear whether the cooling 
in the shocked outer regions would be efficient enough to 
form stars. 

Admittedly, the above arguments do not include com- 
plications such as unusual morphologies of UCHIIs and 
possibly non-spherical ionization front, among others, 
and only provide basic estimates of the gas-removal 
timescales. Observationally, Galactic « 1 Myr old gas- 
free young clusters such as the ONC and the NYC imply 
that the embedded phase is < 1 Myr. The above 
widely-used gas-expulsion model does realize the essen- 
tial dynamical effects on the star cluster. In particu- 
lar, such simplification has practically no effect on the 
remaining cluster after the gas is expelled, e.g., on it's 
re-virialization, which is the focus of this work. 

2.2. Initial configuration, stellar dynamics and 

evolution 

The initial model embedded stellar clusters are Plum- 
mer spheres (Kroupa 2008) which the gas follows as well. 
This is a reasonable approximation since dense interstel- 
lar medium (ISM) filaments appear to have Plummer-like 
sections (Malinen et al. 2012; also see below). The rh{0) 
of embedded clusters are substantially smaller than the 
typical sizes of exposed young clusters as found in the 
semi-analytic calculations by Marks & Kroupa (2012). 
These calculations constrain the birth-density of a large 
number of observed clusters by their observed popula- 
tion of binary stars leading to a remarkable overlap with 
the densities of star-forming molecular clumps (Fig. 6 in 
Marks & Kroupa 2012). Our initial clusters thus follow 
the empirical relation between rh{0) and the embedded 
cluster mass (only the stars) Meci(O) by Marks & Kroupa 
(2012), viz., 

^.o.ois:s;x(^)— , 

which is a rather weak dependence. This independently 

obtained result is in excellent agreement with the ob- 
served results from Herschel (Andre et al. 2011; Fig. 6.3). 
The star formation efficiency is taken to be e = 1/3, i.e., 
Mg{0) = 2Meci(0) (see Sec. 4 for a discussion). 
Notably, Andre et al. (2011) and Malinen et al. (2012) 



refer to the shape and the compactness of ISM filament 
sections. The star forming spherical/spheroidal proto- 
clusters form within the ISM filaments and at their in- 
tersections (Schneider et al. 2012); a part of the filament 
that collapses under self-gravity to form an embedded 
cluster would become spherical/spheroidal. Recent Her- 
schel observations of ridges in molecular clouds and of 
the associated filaments support this (Hennemann et al. 
2012; Schneider et al. 2010; Hih et al. 2011). The com- 
pact sizes (« 0.1 pc) of these ISM fllaments (Andre et 
al. 2011) then dictate the high compactness of the initial 
embedded systems. There is, of course, the competing 
viewpoints whether all the stars form within one cloud 
(the standard scenario) or the final cluster can be formed 
hierarchically (see Sec. 1). Currently, the possibility of 
forming massive stars within the filaments but outside 
any cluster has also been suggested {e.g., Bressert et al. 
2012). So far our computations in this paper are con- 
cerned, we adopt the standard scenario and investigate 
whether the observed kinematics of the yoimg clusters 
like R136 and NYC (to start with) can be reasonably ex- 
plained within such a context. In that case, our chosen 
initial conditions well reflect the sizes and the shapes of 
the observed filaments. 

The IMFs of the clusters are chosen to be canonical 
(Kroupa 2001) with the most massive star following the 
TOmax-Mcci(O) relation of Weidner & Kroupa (2004). To 
that end an "optimized sampling" algorithm (Kroupa et 
al. 2013) is used instead of randomly sampling the IMF. 
All the computed models are fully mass-segregated using 
the method of Baumgardt et al. (2008), as observed in 
young clusters (Littlefair et al. 2003; Chen et al. 2007; 
Portegies Zwart et al. 2010). For computational ease, 
we take all the members of the initial cluster to be sin- 
gle stars. Since the 0-stars in young clusters are typi- 
cally found in binaries (Sana & Evans 2011; Sana et al. 
2012), our initial stellar population, admittedly, does not 
completely represent that of a young cluster. However, 
binarity is unlikely to substantially influence the expan- 
sion of a cluster during gas-expulsion and re-collapse and 
virialization thereafter, which is the focus of this work. 
Although the initial segregated state is detected for sev- 
eral Galactic young clusters, it is yet to be conflrmed 
for young clusters in general (see Portegies Zwart et al. 
2010 for a discussion). Therefore, for comparison pur- 
poses, we compute identical models without any initial 
mass-segregation. We find that mass-segregation does 
not influence the conclusions (Sec. 3.1). 

The dynamical evolution of the model clusters are com- 
puted using the state-of-the-art NBDDY6 code (Aarseth 
2003). In addition to integrating the particle orbits us- 
ing the highly accurate fourth-order Hermite scheme and 



4 



Banerjee and Kroupa 



dealing with the diverging gravitational forces in close 
encounters through regularizations, NBDDY6 also employs 
the well-tested analytical stellar and binary evolution 
recipes of Hurley et al. (2000, 2002), viz., the SSE and 
the BSE schemes. NB0DY6 also includes the time-variable 
Plummer gas potential as described above. 

The initial masses of the computed model embedded 
clusters are chosen to be Moci(O) « IO^Mq which is an 
upper mass limit of R136 (Crowther et al. 2010). Their 
Plummer radii are chosen from Eqn. (2). We also con- 
sider an embedded cluster of Moci(O) = 1.3 x 10^ to 
mimic the gas-removal from a Galactic NYC-like cluster 
(Stolte et al. 2004, 2006; Rochau et al. 2010). The metal- 
licities of the former models are taken to be Z = 0.5Zq, 
as appropriate for the LMC, and that for the NYC model 
it is Z = Zq. Our computed models are summarized in 
Table 1, where those in Kroupa et al. (2001) are also 
added for comparison. For both of our computed mod- 
els, the TgS are similar to the initial crossing times Tcr(O); 
Tg 2 and 1 time(s) Tcr(O) for the R136 and the NYC 
model respectively. Since Tg < Tcr(O) for the Kroupa et 
al. (2001) models, our gas-expulsions are less "prompt" 
owing to higher initial concentrations. The MCLUSTER 
program (Kiipper et al. 2011) is used to set the initial 
configurations. As it is diflScult to model the weak tidal 
field of the LMC, the clusters are evolved in absence of 
any galactic potential. 

3. RESULTS 
3.1. Gas-expulsion from an R136-like cluster 

Fig. 1 shows the evolution of the Lagrange radii, Rf, 
in our computed R136 model clusters (Table 1) with 
Td = 0.6 (top panel) and (bottom) Myr. The = 
case is computed as a test case only to compare with that 
of the more realistic delay = 0.6 Myr. From the be- 
ginning of the gas-expulsion, the cluster expands rapidly 
with timescale Tg. At least 60% of the cluster, by mass, 
then collapses back to a steady size, i.e., regains virial- 
ization in Tvir ~ 1 Myr (beyond Td). The inner regions 
virialize even earlier. Given that the bulk of the R136 
cluster is w 3 Myr old (Andersen et al. 2009), it ought 
to be currently in dynamical equilibrium as a whole even 
if it has undergone a substantial gas- expulsion phase in 
the past. We discuss this further in Sec. 4. The gas-free 
cluster is expanded by a factor of « 3 in terms of its 
half-mass radius, rh, after reverting to dynamical equi- 
librium. Notably, the re-virialization happened in spite 
of the wind mass-loss from, the m,assive stars that was 
operative during these computations. Comparison of the 
two panels in Fig. 1 makes it apparent that a finite Td 
(with Tg unchanged) keeps the form of the cluster's ex- 
pansion and rc-virialization identical — it merely applies 
a time-translation to the overall evolution. 

Fig. 2 shows the i?/-evolution for identical R136 mod- 
els but without initial mass-segregation. It demonstrates 
that the cluster evolutions are throughout identical to 
the previous cases with primordial mass segregation, 
in particular during the expanding phase and the re- 
collapse to virialization. In other words, primordial 
mass-segregation has no effect on the re-virialization of 
a cluster. 

Very recently, Henault-Brunet et al. (2012) mea- 
sured radial/line-of-sight velocities (RV) of single 0-stars 




0.1 : 

0.5 1 1.5 2 2.5 3 

t (Myr) 




0.1 V 



0.5 1 1.5 2 2.5 3 

t (Myr) 

Figure 1. The evolution of the Lagrange radii, Rf, for stellar 
cluster mass fractions / (escaping stars inclusive) for the computed 
R136 models with = 0.6 Myr (top) and r.^ = (bottom). In 
each panel, the curves, from bottom to top, correspond to / = 0.01, 
0.02, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.625, 0.7 and 0.9 respectively. The 
thick solid line is therefore the half mass radius of the cluster. 

within 1 pc < i? < 5 pc projected distance from R136's 
center^ with data obtained from the "VLT-FLAMES 
Tarantula Survey" (VETS; Evans et al. 2011). They 
conclude that the RV dispersion, Vr, of the single O- 
stars within this region is 4 km s^^ < < 5 km s"""^, 
conforming with R136 being in virial equilibrium. 

Fig. 3 shows the evolution of Vr corresponding to the 
R136 computations in Fig. 1 for 1 pc < R < 5 pc 
projection from the density center and stellar masses 
M > 16Af0, i.e., those of the 0-stars. It demonstrates 
that, after re-virialization, the 0-stars indeed have Vj. 
remarkably similar to that meastired by Henault-Brunet 
et al. (2012). The initial large fiuctuations in Vr in the 
upper panel of Fig. 3 (r^ = 0.6 Myr) is due to the initial 
mass-segregated condition which results in only a few O- 
stars within 1 pc < i? < 5 pc initially. The measured RV 
data is therefore consistent with R136 being through a 
gas-expulsion phase and then in a re-virialized state now, 

^ Strictly, the distances were measured from the most massive 
cluster member star R136al. The exact location of the R136's true 
density center being unknown, we consider this star be expectedly 
very close to the cluster's density center. 
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Figure 2. The evolution of Lagrange radii, Rr, for computed 
models without any initial mass-segregation but otherwise identical 
to the R136 models in Fig. 1. The legends are the same as in Fig. 1. 

unlike claimed by Henault-Brunet et al. (2012). 

3.2. The case of the NGC 3603 Young Cluster 

Fig. 4 (top) shows the Lagrangian radii of our com- 
puted NYC cluster (Table 1). The mass of NYC (Stolte 
et al. 2004, 2006; Rochau ct al. 2010) is substantially 
smaller than that of R136, implying a longer free-fall and 
hence re-virialization time. The model is computed with 
Td — 0, which is rather unrealistic, to see the effect of 
the earliest possible gas removal. Even then the cluster 
is still expanding in its outer parts and just beginning to 
fall back in the inner regions (< -R0.3) at its very young 
age of « 1 Myr (Stolte et al. 2004), as can be seen in 
Fig. 4. In other words, it is practically impossible for 
NYC to be in dynamical equilibrium at its present age 
had it undergone a substantial gas-expulsion phase ear- 
lier. The computed cluster re-virializes in Tvir > 2 Myr 
(< 40% of it; c.f. Fig. 4) as opposed to our R136 models 
(Sec. 3.1) which take « 1 Myr for the same (> 60% of 
them; c.f. Fig. 1, Table 1). 

From proper motion measurements within the central 
« 1 pc projection of NYC, Rochau et al. (2010) found no 
notable differences among the 1-dimensional velocity dis- 
persions (Vid) of 4 mass (magnitude) bins over 1.7— QMq. 
This indicates that NYC is far from energy equipartition 
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Figure 3. The evolution of the radial velocity (RV) dispersion, Vr, 
of the O-stars (M > 16Mq), within the projected distances 1 pc < 
-R < 5 pc from the cluster center, for the computed R136 models. 
The top and the bottom panels correspond to the computations 
with Tcj = 0.6 and Myr respectively. 

as it should be if the system is out of dynamical equi- 
librium (the converse conclusion that non-equipartition 
implies non-equilibrium, as used by Rochau et al. 2010, is 
however not necessarily true). This is explicitly demon- 
strated in Fig. 4 (bottom), where the initial separation 
in Vid is practically vanished at t = 1 Myr due to violent 
relaxation driven by the gas-expulsion. 

The computed Via at t = 1 Myr is somewhat less than 
the measured^ Vid ~ 4.5 ± 0.8 km s""'^ and the observed 
system must currently be super-virial, unlike the com- 
puted sub-virial state at t = I Myr, as the inferred dy- 
namical mass^ exceeds the photometric mass (Rochau et 
al. 2010). These can be easily accoimted for by a plau- 
sible 0.6 < Td < 0.8 Myr delay in the gas removal {c.f. 
Fig. 4) which would then result in agreement with the 
observed Vid and the super-virial state of NYC. Further- 
more, NYC's mass is uncertain by a factor w 2 (Rochau 
et al. 2010). 

The is obtained from proper motion measurements and is 
therefore "binary-corrected" . 

^ The dynamical 1.7 X W^Mq, doesn't represent an upper limit 
as quoted by Rochau et al. (2010), but is only an individual esti- 
mate. 
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Figure 4. Top; The evolution of the Lagrange radii for the com- 
puted NGC 3603 Young Cluster model with = 0. The leg- 
ends arc the same as in Fig. 1. Bottom: The variation of the 
1-dimensional (transverse) velocity dispersion, Vid, within R < 1 
pc projection, for 4 mass-bins chosen within the same mass rajige 
as in Rochau et al. (2010), for the above model. Note that initially 
Vid differs between the mass bins consistent with energy equipar- 
tition. This separation of Via values vanishes as a result of violent 
relaxation driven by gas-expulsion. Note that the t-axis of the 
bottom panel is plotted in the logarithmic scale to highlight this. 

Fig. 5 replots Fig. 4 with a « 0.6 Myr delay which, 
in turn, corresponds to the resuhs of a computation with 
this delay in gas-expulsion (see Sec. 3.1). This makes the 
cluster super- virial at its « 1 Myr current age. The corre- 
sponding Vid is still somewhat smaller than the observed 
value but is within 3(t limits. Notably, NYC might, in 
fact, be somewhat younger (c./. Fig. 4 of Stolte et al. 
2004) and Vid also depends on the exact initial mass and 
size which are subject to uncertainties. Moreover, NYC's 
Vid and age are susceptible to its distance uncertainties. 
Notably, Cottaar et al. (2012) also find the Westerlund I 
cluster sub-virial. 

4. CONCLUSIONS AND OUTLOOK 

The key message from our above computations and 
analyses is a reminder that an observed dynamical equi- 
librium state of a very young stellar cluster does not 
necessarily dictate that the cluster has not undergone a 
gas-expulsion phase. The non-occurrence of gas-removal 



Q I I I I I I I 

1 1.5 2 2.5 3 3.5 

t (Myr) 

Figure 5. Evolutions of the Lagrange radii (top) and 1- 
dimensional velocity dispersion (bottom) of the same computed 
NYC models as in Fig. 4 where an initial gas expulsion time-delay 
of = 0.6 Myr is applied while plotting. Hence, the plots are 
essentially the same as in Fig. 4 except that they begin at t = 0.6 
Myr and correspond to those from a computation with = 0.6 
Myr (see Sec. 3.1). They show that NYC would be super-virial at 
its current Ri 1 Myr age, as observed. 

has been highlighted, particularly by putting forward the 
case of R136, by recent authors, e.g., Henault-Brunet 
et al. (2012). This is why we focus on R136 in this 
study. We also demonstrate, choosing the example of 
NGC 3606, that lower-mass clusters take much longer to 
re-virialize (c./. Table 1) so that it is as well appropri- 
ate to find non-equilibrium signatures in them. This has 
been shown to be true also for the ONC by Kroupa et al. 
(2001) (c./. Table 1). The above must be remembered 
while interpreting the kinematics of any young cluster. 

Notably, R136 exhibits an age-spread; its massive stel- 
lar population can be younger, viz., < 2 Myr old (Massey 
& Hunter 1998). If the gas-expulsion is assumed to com- 
mence only after the formation of the most massive stars, 
a currently equilibrium state is still plausible given the 
Tvir ~ 1 Myr re-virialization time. It should, however, 
be remembered that age-estimates of massive stars are 
highly uncertain. There isn't any concrete evidence that 
invalidates the scenario involving the formation of R136's 
entire stellar population at once « 3 Myr (Andersen et 
al. 2009) ago and the most massive single stars forming 
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later via dynamically induced binary mergers (Banerjee 
et al. 2012b) or appearing younger due to mass transfer 
in close 0-star binaries. 

Although the present study does not cover a general 
range of the parameters r^, Tg, r/j(0) and e, these com- 
putations can still provide a fair idea of the effects of the 
variation of these parameters over their plausible ranges. 
By applying time-shifts to the point of expansion of the 
clusters in Figs. 1 & 2, it can be concluded that < 2 
Myr required to find R136 virialized at t = 3 Myr, since 
Tvir ~ 1 Myr. Although an unambiguous ra requires 
radiation hydrodynamic calculations, this upper limit is 
perhaps too large a value for in the light of the discus- 
sions in Sec. 2.1. The currently chosen values of TgS are 
upper limits based on the typical H Il-gas sound-speed of 
w 10 km s~^ (Sec. 2.1). As discussed in Sec. 2.1, the out- 
flow speed can initially be significantly higher than the 
sound-speed, driven mainly by radiation pressure, result- 
ing in smaller TgS than those assumed. This would make 
the gas-expulsion more prompt, keeping Tvir practically 
unchanged. 

As for the SFE, the chosen e w 1/3 is in the range 
of the observationally estimated SFEs by Lada & Lada 
(2003) (see their Table 2). Recent theoretical modelling 
of star formation with high-resolution resistive magneto- 
hydrodynamics (Machida & Matsumoto 2012) also s\ig- 
gest e « 33%. In any case, a larger e would also result 
in a bound system and shorter Tvir as the re-virializing 
mass would be larger. As for r?i(0)s, their increasing val- 
ues would increase TvirS, so our conclusions depend on our 
adopted initial high compactness. However, this condi- 
tion is plausible since such compact star-forming envi- 
ronments (a factor of w 10 smaller than typical young 
clusters, i.e., rh{0) « 0.3 — 0.4 pc) are inferred obser- 
vationally by Andre et al. (2011) and as well from the 
semi-analytic study by Marks & Kroupa (2012) indepen- 
dently (the re-virialization time depends on the above 
typical compactness but, of course, not on the specific 
Eqn. (2)). In other words, plausible variations of the 
model parameters are unlikely to alter our above primary 
conclusions since none imply a substantial lengthening of 

Td, Tg or Tvir. 

An immediate improvement over the present study 
would be to introduce an appropriate binary popula- 
tion as in Banerjee et al. (2012a, b). In the light of the 
present and upcoming observations of Galactic and ex- 
tragalactic yoking clusters, an important pending task is 
to systematically and quantitatively study the effect of 
gas-expulsion over their entire mass-range with varying 
parameters such as r^, Tg, rh(0) and e including real- 
istic stellar populations. Such a model-bank would be 
extremely valuable to interpret the ever-enriching kine- 
matic data of young stellar clusters. 
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